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Summary 

A typical problem in optimal design theory is finding an experimental design 
that is optimal with respect to some criteria in a class of designs. The most 
popular criteria include the A- and D-criteria. Regular graph designs occur in 
many optimality results and, if the number of blocks is large enough, they are 
A- and ID-optimal. We present the results of an exact computer search for the 
best regular graph designs in large systems for up to 20 points, k < r < 10 and 
r{k — 1) — (u — 1) [r{k — l)/{v — 1)J < 9. 
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1 Introduction 


Suppose we are to design the following statistical experiment: there are v treatments to 
be compared on a number of experimental units that can be partitioned into b blocks 
of size k with k < v. Typically, the blocks might differ systematically but all units in 
a block are assumed to be alike. For hxed u, b and fc, how should the treatments be 
allocated to the units to get as much information as possible from the available data? 
The latter often means the estimate of the unknown parameters with the least possible 
variance. If there are several parameters, this is a multidimensional problem and the 
design of the experiment can be ‘good’ in different ways. 


We will give the basic dehn i tions in the following and the reader is referred to the 
paper by iBailev fc Cameron! (1200911 . where a good overview and more details on the 
application of combinatorics can be found. Formally, a block design d is an assignment 
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of V treatments or varieties to a set of experimental units that have been partitioned into 
b blocks of size k. The statistical model to analyse the data obtained from d is assumed 
to be the linear model. For each yield of unit u in block j, the model is 

Yu)j = Tf(^coj) + fdj + 

where (3j is the effect of the block j, Tj is the effect of treatment i and e^ij are uncorrelated 
random errors with expectation 0 and variance The function f{ujj) in the above 
equation specihes which treatment is allocated to unit u in block j. 

The design is said to be connected if all pairwise differences Tj — tj are estimable. A 
design is called complete if every treatment occurs in every block; otherwise it is called 
incomplete. If every treatment occurs at most once per block, the design is called binary 
and we will always assume that all designs we consider are incomplete, binary and 
connected. The replication ri of a treatment i is the total number of units that have 
been assigned treatment i. If the replications are all equal to a constant r, the design is 
called equireplicate. A well-known example for block designs are the balanced incomplete 
block designs (in short BIBDs), these are binary equireplicate incomplete block designs 
with replication bk/v, where k < v and any pair of treatments is contained in exactly A 
blocks for some A > 0. BIBDs with these parameters are also called 2-{v, k, A)-designs. 

Let Nd be the v xb treatment-block incidence matrix of the design d, i.e. the ij-entry of 
Nd is the number of units in block j that have been assigned treatment i. The ij-entry 
of the product NdNj will be denoted by A^- and is called the concurrenee of treatments 
i and j. Note that for binary designs Xu = r* and Xij is the number of blocks containing 
both treatments i and j. The information matrix for the estimation of the treatment 
effects is 

Cd = diag(ri,..., r^) - -j^NdNj. 

Since Cd has row sums 0, the all-1-vector is an eigenvector with eigenvalue 0. All other 
eigenvalues are called the non-trivial eigenvalues. Note that if the design is connected, 
the rank of the information matrix is u — 1, the eigenvalue 0 has multiplicity 1 and all 
non-trivial eigenvalues are strictly positive. We will use another matrix associated with 
the design: The Laplacian matrix of a binary connected design is 

Ld = k diag(ri,..., r„) - NdNj = kCd. 

Let pi{d) > P 2 {d) > ... > py_i{d) > Pv{d) denote the eigenvalues of Ld. Note that if 
the design is connected, then p^_i{d) > Pv{d) = 0. The eigenvalues pi((i),..., 
are called the non-trivial Laplaeian eigenvalues. 

The average variance of the set of the best linear unbiased estimators of the pairwise 
differences of the treatment effects is proportional to the reciprocal of the harmonic 
mean of the non-trivial eigenvalues of Cd. The volume of the conhdence ellipsoid for 
the estimate of (ti, ... ,Tv) is proportional to the square root of the reciprocal of the 
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product of the non-trivial eigenvalues of Cd (for example Shah &: Sinha fllQSg lh These 
facts give us different criteria with regards to which a design can be considered ‘good’: 
A design maximizing the harmonic mean of the non-trivial eigenvalues of Cd is called 
A-optimal and a design is called D-optimal, if it maximizes the product of the non-trivial 
eigenvalues of Cd- There are many more optimality criteria; another popular example is 
the i?-criterion which is the maximization of the smallest non-trivial eigen value of Cd and 
i s equ ivalent with minimizing the largest variance of the estimators (see 
(Il989h b 


Shah fc Sinha 


There is no general answer to the question of which design is to b e cho sen for given 
{v,b,k), but there are several partial results. For example, iKieferl fll975l) proved that 
BIBDs are optimal with regards to a wide range of criteria, in particular the A- and 
D-criteria. But it is not clear which designs to choose if no BIBD exists. One class of 
designs that have been suggested to be good candidates are the regular graph designs 
(in short RGBs). RGBs are equireplicate binary designs that are ‘close’ to BIBBs in the 
sense that any pair of points occurs in either A or A -|- 1 blocks for some integer A > 0. 
The name hints towards a connection with regular graphs which we will explain in the 
hrst section. The Laplacian matrix of an RGB d with v points, replication r and block 
size k can be written as 


Ld = {r{k - 1) -h A}G -Td- XJv, 


where G is the v x v identity matrix, Jy denotes the v x v all-1-matrix and Td is a 
symmetric n x n (0, l)-matrix with O’s on the diagonal and exactly r{k — 1) — X{v — 1) 
number of I’s in each row and each column. Note that if A = 0, the fact that the design 
is connected is equivalent with the matrix Td being irreducible (i.e. the matrix can not 
be transformed into a block upper-triangular matrix by row or column permutations). 


John fc Mitchelll ( 1977 1 conjectured that if an incomplete block design is B-optimal (or 


A-optimal or F'-optimal), then it is an RGB (if any RGBs exist). In the same paper, they 
provided a list of the best RGBs with n < 12, r < 10 and v < b which they found using 
n umeri c al me thods. RGBs occur as optimal designs in many results, see for example 


Ghena fll98ll) and are often candidates to be used when designing an experiment in 
small cases. Therefore, a list of best RGBs will be very useful in practical design of 
experiments. 


Howev er, John and Mitchell’s conjecture is not true in gene ral (see for ex ample iBailev 
( 2007 1). but it holds if the number of blocks is large enough ( Gheng ( 19921) 1. Numerical 
methods as used by I John fc Mitchelll (119771) will not extend to these large systems, but 
we can use ex act methods capital i zing o n the power of symbolic algorithms provided by 
Mathematica. I Wolfram Research! (1201211 . 


Suppose, d is a 2-(n, fc, A)-design on v points and block size k with b blocks and Laplacian 
matrix L. Then for y eN the matrix L[y] = Ld + yL^ is the Laplacian matrix of a design 
on V points, replication r + yX{v — l)/{k — 1) and b + yb blocks of size k. We can write 
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L[y] = Ld + yLj with A = [r{k — l)/{v — 1)J and 5 = r{k — 1) — A(n — 1) as 


L[y\ = {{r{k-l) + \]I^-Td-\Jv + yK'vIv-Jv) 

= {6+ vy\ + vX)I^-{y\ +X)J^-Td 
and with x = \ + y\ we have 

L[x\ = ((5 + vx)I^ — Td — xJ^. 

To identify the best RGD for large y, we take a similar approach to John and Mitchell: 
we order all possible matrices T according to the performance of the corresponding L[x] 
on the A- and D-criteria for different values of x. In particular, we hnd values Xq,Xq 
such that these orders do not change for all x > Xq,Xq . To hnd the best RGBs with 
regards to the A- and B-criteria for given r and k, all one has to do is to hnd the hrst 
matrices in the orders that produce a design. Since we are only searching among RGBs, 
we will say that these designs are the A- or B-best RGBs. 

We would like t o high light some key diherences between our approach and that of 
John fc Mitchelll (1977). Firstly, we have only considered A- and B-optimality since 
there is a list of binary connected B-optimal designs u p to 15 points without the re¬ 
striction to RGBs due to a construction by iMorgaiJ (120071) which can be found on 
www.designth eory.org. S econd ly, to obtain a list of all the matrices T we use the pro¬ 
gram genreg flMeringerl fjl999l) h It is due to the restrictions of this program that we 
can only consider irreducib l e mat rices with row sum 6 < 9. The fact that all optimal 
designs in iJohn fc Mitchelll (119771 1 with S > 2 correspond to irreducible matrices makes 
it a reasonable pay-off to be able to handle cases with a large number of matrices. We 
like to note however that RGBs corresponding to irreducible matrices might not per¬ 
form better than the ones corresponding to reducible matrices, but we failed to hnd an 
example. Also note that there is only one irreducible matrix for J = 2orJ = n — 1 and 
we therefore exclude these cases. And hnally, as mentioned above, we are using exact 
methods and therefore the orderings of the matrices are valid for all x > xa,xd- 

Our results conhrm that all the designs listed in John &: Mitchelll ( 19771) are indeed 
the best RGBs and moreover, that they are the best RGBs for big y (due to exact 
computation). We also found A- and B -best RGBs for v = 11, k = r = 3 and k = r = 8] 
these are two cases that are missing in iJohn fc Mitchelll (119771) . We extended the list of 
the best RGBs to n < 14 and k < r < 10, S < 9 for all admissible block sizes and up to 
n < 20 and k = 2. 


A list of all found RGBs not listed in iJohn fc Mitchelll (119771 ) can be found in the ap¬ 
pendix and a full list of all found designs in xm l format is inc luded in the electronic 
appendix. For more information on xml hies see ISoicherl (120061) . The search produced 
examples supporting some open conjectures, which we discuss in the hnal section. More¬ 
over, we have found an example with n = 14 and y = 0 where the A-best RGB is not 
B-best and vice versa. 
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2 Preliminaries 

We will denote the set of all symmetric irreducible u x u (0, l)-matrices with O’s on the 
diagonal and row and column sum S by Ai{v,S). 

A graph Q is given by a set of vertices V{Q) and a set E{Q) of two-element subsets of 
V{Q), called edges. We allow that two vertices are joined by multiple edges but assume 
that there are no loops, i.e. there is no edge of the form {u,u} for any u G V{Q). If 
Q contains no multiple edges, we say that Q is simple. The number of edges / G E{Q) 
with M G / for a vertex u G V{Q) is called the degree of u. If the degrees of all vertices 
of Q are all equal to a constant h, then Q is called h-regular. A graph can be expressed 
as a matrix as follows: the adjacency matrix A(Q) is the symmetric |U(^)| x |U(^)| 
matrix whose ij-entry is the number of edges joining vertices i and j. Note that if Q is 
(5-regular, then any row and any column of A(Q) sums to S. If any vertex of Q can be 
reached from any other vertex by going along edges, then Q is called connected. In this 
case, A(Q) is irreducible and if ^ is simple A{Q) G M(v,S). 

Every design d with parameters (u, b, k) gives rise to a graph in a natural way by taking 
the treatments as vertices and joining any two distinct vertices i and j by Xij edges, where 
Xij is the concurrence of treatments i and j. This graph is called concurrence graph and 
we will denote it by Qd- Note that Qd contains no loops. The design is connected if and 
only if the graph Qd is connected. If the design is binary and equireplicate with replication 
r, then Qd is the r(A; —l)-regular graph with adjacency matrix A(Q) = NdNj — r{k — l)I^. 
The RGBs owe their name to the fact that their adjacency graph is obtained from a 
simple regular graph Q by adding A edges between any two pairs of vertices. Here, 
A = \r{k — l)/{v — 1)J, the graph Q has degree 5 = r{k — 1) — X{v — 1), its adjacency 
matrix is a (0, l)-matrix Td and the Laplacian matrix of the RGB can be written as 

Ld = {r{k - 1) + A}/„ - Td - XJy. 

That means, the Laplacian matrix of any RGB is fully determined by the parameters v, 
k, r and the matrix T^. In particular for A = 0, the matrix Td is the adjacency matrix 
of a connected graph and is therefore irreducible and Td G Ai{v, (5). 
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Let d be an RGD with Laplacian matrix with v points and block size k with b 
blocks and replication r. Fnrther, let L^he the Laplacian matrix of a 2-{v, k, A)-design 
with b blocks, where A is minimal snch that a 2-(n,/c, A)-design exists. We can write 
L[y] = Ld + yL^ with A = \r{k — l)/(n — 1)J, 5 = r{k — 1) — \{y — 1) and x = A + |/A as 

L[x] = (5 + nx)4 - Trf - xJy. (1) 

Let 'il^i > > ■ ■ ■ > 'il^v denote the eigenvalnes of 61^ — Td. Then, ■0^ = 0 since Td 

has row snm 6. The non-trivial eigenvalnes of L[x] are vx + 'ijji,... ,vx + Let i/’t 

denote the vector (-01,... The prodnct of all non-trivial eigenvalnes of L[x] can 

be written as a polynomial in x: 

V—1 V—1 

x) = JJ(nx + tl^i) = '^{vxy~^~^Sj{'il)T), (2) 

i=l j=0 


where Sj denotes the jth elementary symmetric polynomial, that is S'o = 1 and for 
z = (zi,... ,Zy-i) e 

5'j(z) = j = 

The harmonic mean of all non-trivial eigenvalnes of L[x] is 

n — 1 (v — l)iA('i/? 2 ., x) 


A{xPt,x) = 




With the relation 


V—1 




i=i 


(see for example iBeckenbach &: BellmanI fjl965l) . p. 34) we have 

V—1 V—1 V—1 v—2 

JJ {vx + yy = 

i=l,i^l 1=1 j=0 

v-2 

j=0 

Let Dx{%l>rp^x) = ^iA('i/?y,x) denote the derivative of D('ipj,,x) in x, then 


A{tJjrp^ x) = v{v — 1) 


Dy'ijjT,x) 


(3) 


Note that in particnlar the Laplacian matrix Ld of the RGD is Ld = L[A] and D{'ij)rp^ A) 
and A{'il)rp^ A) are the prodnct and harmonic mean of its non-trivial eigenvalnes respec¬ 
tively. 
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3 Computation method 


For the matrix L[x\ to be the Laplacian matrix of an existing design, only some valnes for 
X will be admissible. Bnt comparing designs with Laplacian matrix L[x\ is now rednced 
to comparing the valnes D{xl)j,,x) and A^'iprp^x) among the matrices T G Ai{v,S). 

The direct approach of compnting the eigenvalnes of L[x] and compnting A{xl^rp^x) and 
D{'if;rp^x) as their harmonic mean and prodnct leads to long compntation times if exact 
methods are used. It is possible however to compute A{il)rp^x) and Di^rp^x) more 
efficiently: Let pL = (pi, p 2 , • • •, P^-i) denote the vector of the non-trivial eigenvalues of 
L[x\. The characteristic polynomial x{^) of L[x] can then be written as 


X{x) = 


V 

j=0 


-iySj{pL)x 


v-j 


With equation[2l we can compute x) with exact methods in Mathematica. IWolfram Research 

(120121) as the last non-zero coefficient of x(a;). Since we can express A{xl)p, x) as in equa- 
tion[31 it is enough to compute D{'il)p,x). 


To find the best RGBs on v points, we first generate all matrices in J\A{v, 5). We do this 
by generating the adjacenc y lists of all connected 5-regular graphs on v vertices with the 
program genreg (IMeringeii (119991) 1 for 5 < 9 from which we compute the matrices L[x] 
as given in equation [H As explained above, it is enough to order all matrices in A4{v, 6) 
corresponding to the value A) and the value A^y^p, A) with A = Lr(fc —l)/(n —1)J. 


We will explain this in more detail in case of the value A(V’r 5 fh® same procedure is 
applied for finding the orders with regards to Dl'if^p, A): To find the best RGB, we want 
to order the matrices, say A4(v, 5) = {Ti,..., Tm} such that 


A^'ij^p.jx) > A{'if)p.^^,x) for all a; > A and z = 1,..., M — 1 
or, equivalently such that the polynomial 


P{Ti,Ti+i,x) = D{'ipp.,x)D^{'ipp.^^,x) - D{'if;p.^^,x)D^{'if;p.,x) 
has no roots bigger than A, that is if P{x) = 0, then x < X. 

Here, A still depends on the existence of RGBs corresponding to a matrix in Ai{v,6). 
Instead of computing the exact values for A, we want to find the smallest integer x^ 
(not depending on the existence of designs) and the order of the matrices such that the 
following two conditions are satisfied, in which case we say that the order stabilizes for 

rj.A 

Xq . 

1. There exists an x < x^ such that 

P(Ti, Tj+i, x) < 0 for z = 1,..., M - 1; 


and 
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2. for all X > Xq, 


P{Ti, Ti+i,x) > 0 for z = 1,..., M - 1. 


That means, if the order of the matrices stabilizes for x^, we have found the order of 
the matrices according to A{'iprp^ A) for any admissible A > Xq. We will say that Ti has 
rank i in this order. 


We denote by x^ the equivalent of Xq for the order according to for T G 

M{v,5). 


We obtain these values by first guessing a value for Xg and Xq and then verifying the 
above two properties for all matrices in A4(v, 6). There exist exact alg orithms for the 
search of roots of a polynomial in Mathematica, IWolfram Researchl (120121) . The following 
table shows the values of Xq and Xq for different v and 6. 


V 

6 

Xq 

rfD 

Xq 

5 

4 

0 

0 

6 

3 

0 

0 


4 

0 

0 

7 

4 

0 

0 


6 

0 

0 

8 

3 

0 

0 


4 

0 

0 


5 

0 

0 


6 

0 

0 


V 

5 

Xq 

Xq 

9 

4 

1 

1 


6 

0 

0 

10 

3 

1 

1 


4 

1 

1 


5 

1 

1 


6 

1 

1 


7 

0 

0 


8 

0 

0 

11 

4 

1 

1 


V 

5 

^A 

Xq 

rfD 

Xq 

11 

6 

1 

1 


8 

0 

0 

12 

3 

1 

1 


4 

2 

2 


5 

3 

3 


6 

3 

2 


7 

2 

2 


8 

1 

1 


9 

0 

0 


V 

5 

r^A 

Xq 

rfD 

Xq 

13 

4 

3 

3 


6 

4 

3 


8 

5 

3 

14 

3 

1 

1 


4 

4 

3 


5 

6 

5 

15 

4 

5 

4 

16 

3 

1 

1 

18 

3 

2 

1 

20 

3 

2 

2 


To find the best design with replication r and block size fc, we search for the first 
matrix among the ordered m atrices t h at giv es rise to a block design. To do this we use 
the GAP package DESIGN (Soicher f 20061) b Note that in some cases, due to either 
the large number of graphs or long computation times with GAP, we could not search 
among all matrices in Ai{v,6) and this might be a reason why we could not find any 
designs for some choices of r and k. 


3.1 Table of A- and D-hest RGDs in large systems 


There follows a table of the A- and D-best RGBs (for y > max{Cf (x^ — A)/A}). We list 
6 = r{k — 1) — A(u — 1) and the smallest A such that a 2-(u, k, A)-design exists (found 
with GAP). 


Most designs can be found in either IClatworthvI (119731) or iJohn et al.l (119721) . To follow 
the convention as in iJohn fc Mitchell! (119771) . we write in the reference column P.XY 
for a design in iGlatworthvi (119731) with reference number XY and G.XY for a design in 
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John et al.l (119721) with reference number XY. If the design is not in either catalogues 


but can be found in I John fc Michelll (119761) we give the reference number as JM.XY and 
if it is cyclic we give the initial blocks. If the design is the complement of a design in 
one of the catalogue, we add an (C) to the reference number. All other designs can be 
found in the appendix. However, so me of these designs are possible to construct with 
known methods, such as the ones in l.lohni (1196711 . 


V 

k 

r 

A 


Reference 

V 

k 

r 

A 


Reference 

6 

2 

3+5y 

1 

3 

RSR6, C.A2 

9 

3 

10+4y 

1 

4 

JM.IO 

6 

2 

4+5y 

1 

4 

RR18, C.A3 

9 

4 

4+8y 

3 

4 

C.B12 

6 

2 

8+5y 

1 

3 

RR24 

9 

5 

5+lOy 

5 

4 

C.B12(C) 

6 

2 

9+5y 

1 

4 

RR27 

9 

6 

6+8y 

5 

6 

RSR65 

6 

3 

4+5y 

2 

3 

(1,2,4)(1,3,5) 

10 

2 

3+9y 

1 

3 

RT2 

6 

3 

7+5y 

2 

4 

P.R46 

10 

2 

4+9y 

1 

4 

C.A16 

6 

3 

9+5y 

2 

3 

RR52 

10 

2 

5+9y 

1 

5 

RSRll, C.A17 

6 

4 

6+lOy 

6 

3 

RSR35 

10 

2 

6+9y 

1 

6 

JM.12 

6 

4 

8+lOy 

6 

4 

RR96 

10 

2 

7+9y 

1 

7 

JM.13 

7 

2 

4+6y 

1 

4 

#1 

10 

2 

8+9y 

1 

8 

RR36, C.A20 

7 

2 

10+6y 

1 

4 

JM.l 

10 

3 

3+9y 

2 

6 

JM.14 

7 

5 

10+15y 

10 

4 

#2 

10 

3 

6+9y 

2 

3 

RT12 

8 

2 

3+7y 

1 

3 

C.A7 

10 

4 

4+6y 

2 

3 

RT33 

8 

2 

4+7y 

1 

4 

RSR9, C.A8 

10 

4 

8+6y 

2 

6 

JM.15 

8 

2 

5+7y 

1 

5 

JM.3 

10 

4 

10+6y 

2 

3 

RT37 

8 

2 

6+7y 

1 

6 

RR29, C.AIO 

10 

5 

6+9y 

4 

6 

JM.16 

8 

2 

10+7y 

1 

3 

(1,2)(1,2)(1,3) 

10 

5 

8+9y 

4 

5 

JM.17 






(1,4)(1,5)(1.5) 

10 

5 

10+9y 

4 

6 

JM.19 

8 

3 

3+21y 

6 

6 

RR54, C.B3 

10 

6 

6+9y 

5 

3 

RT60 

8 

3 

6+21y 

6 

5 

JM.4 

10 

7 

7+21y 

14 

6 

JM.14(C) 

8 

3 

9+21y 

6 

4 

RR58 

11 

2 

4+lOy 

1 

4 

C.22 

8 

4 

4+7y 

3 

5 

C.B6 

11 

2 

6+lOy 

1 

6 

#3 

8 

4 

6+7y 

3 

4 

RSR38 

11 

2 

8+lOy 

1 

8 

JM.20 

8 

4 

8+7y 

3 

3 

(1,2,3,5)(1,2,3,6) 

11 

3 

3+15y 

3 

6 

#4 

8 

4 

9+7y 

3 

6 

RRlOl 

11 

3 

9+15y 

3 

8 

JM.21 

8 

5 

5+35y 

20 

6 

RR314, C.B3(C) 

11 

4 

8+20y 

6 

4 

C.B25(C) 

8 

5 

10+35y 

20 

5 

JM.4(C) 

11 

8 

8+40y 

28 

6 

#4(C) 

8 

6 

9+21y 

15 

3 

C.A7(C) 

12 

2 

3+lly 

1 

3 

JM.22 

9 

2 

4+8y 

1 

4 

C.All 

12 

2 

4+lly 

1 

4 

C.A27 

9 

2 

6+8y 

1 

6 

RR34, C.A13 

12 

2 

5+lly 

1 

5 

#5 

9 

3 

3+4y 

1 

6 

RSR23 

12 

2 

6+lly 

1 

6 

RSR13,C.A29 

9 

3 

6+4y 

1 

4 

JM.8 

12 

2 

7+lly 

1 

7 

JM.24 

9 

3 

7+4y 

1 

6 

RR62 
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V 

k 

r 

A 

5 

Reference 

12 

2 

8+lly 

1 

8 

RR38, C.A31 

12 

2 

9+lly 

1 

9 

P.R39, C.A32 

12 

3 

3+lly 

2 

6 

#6 

12 

3 

4+lly 

2 

8 

P.SR26 

12 

3 

7+lly 

2 

3 

JM.25 

12 

3 

9+lly 

2 

7 

JM.27 

12 

3 

10+lly 

2 

9 

P.R78, C.B34 

12 

4 

3+lly 

3 

9 

P.SR41 

12 

4 

5+lly 

3 

4 

C.B37 

12 

4 

6+lly 

3 

7 

JM.28 

12 

5 

5+55y 

20 

9 

P.R145, C.B43 

12 

6 

10+lly 

5 

6 

P.SR71 

12 

7 

7+77y 

42 

9 

P.R176, 

C.B43(C) 

12 

9 

9 +33y 

24 

6 

#6(C) 

13 

2 

4+12y 

1 

4 

C.A35 

13 

2 

6+12y 

1 

6 

#r 

13 

2 

8+12y 

1 

8 

#8 

13 

3 

9+6y 

1 

6 

#9 

13 

5 

5+15y 

5 

8 

C.B58 

14 

2 

3+13y 

1 

3 

#10 

14 
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4+13y 

1 

4 

#11 

14 

2 

5+13y 

1 

5 

# 12 

14 

3 

9+39y 

6 

5 

# 13 

14 

3 

15+39y 

6 

4 

#14 

14 

6 

6 +39y 

15 

4 

C.B72 

15 

2 

4+14y 

1 

4 

#15 

15 

3 

9+28y 

1 

4 

C.B77 

16 

2 

3+15y 

1 

3 

#16 

16 

3 

9+13y 

2 

3 

#17 

16 

4 

6+5y 

1 

3 

P.R118 

18 

2 

3+17y 

1 

3 

#18 

18 

3 

10+17y 

2 

3 

#19 

18 

5 

5 +85y 

20 

3 

#20 

20 

2 

3 +19y 

1 

3 

#21 
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3.2 Summary and remarks 

Comparison with the results of John and Mitchell 

For all V < 12 and 5 < 5, the best RGBs we found for y = D are isomorphic to the 
designs presented by John and Mitchell except for the cases n = ll,/c = r = 3 and 
k = r = 8. All other designs listed by John and Mitchell for 6 > 2 correspond to the 
same matrix in M{v,6) as the best RGBs we found. This is particularly remarkable 
since we made the additional restriction on all matrices to be irreducible. In the cases 
V = 11, k = r = 3 and k = r = 8 we found an A- and B-best RGB that John and 
Mitchell failed to hnd. The A- and B-best RGBs correspond in both cases to the same 
matrix in Af (11, 6). 


Order stabilization and ranks of the matrices in Ai{v,6) 

In all cases we found that if T G A4{v,6) has rank i in the order corresponding to 
Ai^jprp, Xq), then it has rank i in the order corresponding to D{'if)rp, Xq). The same is not 
true in general for the orders for x > 0; in particular it is not true for any case where 
V > 13. 

The matrix with rank 1 in both of the orders for x = 0 remains A- and B- best with 
growing x, except for the case x = 14, 5 = 5. In this case, the same is true for the order 
of D{'ij)rp,x), but not of A{'ipj,,x) (see also Remark [1]). 

And hnally, we have Xq,Xq <5 + 1; the worst case is xj^ = 6 and x^ = 5 for x = 14, 
5 = 5. 


A- and B-best RGDs 

In most cases the best design was found with a matrix among the best two matrices and 
except for only a few cases, the matrix was among the best ten. 

Except for the case x = 14, 5 = 5 the best RGBs for both A{'ij)rp,0) and B('0^,O) 
correspond to the same matrix in A4(x, 5). Note that since the ranks of the matrices in 
the orders corresponding to A{'iPj,,Xq) and D{'iprp,XQ) are the same, the best RGBs for 
both A{xI)j,,Xq) and D{'iPj,,Xq) correspond to the same matrix in M.{v,5) in all cases. 

Remark 1. The case x = 14, r = 5, /c = 2: The best matrix in A4(14, 5) in the order 
for B('0/p,O) (which is the same for all x > 0) is only second best for ^(' 02 ^, 0 ), but 
becomes best for A{'il)j<,x) where x > 1. The best matrix for A(i/>^,0) is second best 
for A{'il)rp,x) where x > 1. This means in particular that the A-best RGB for x = 14, 
r = 5 and k = 2 for ?/ = 0 is not B-best and vice versa. 
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The y4-best RGD for j/ = 0 is the following design: 


1 2 
2 10 
5 7 

8 13 

13 14 


1 3 

3 7 

5 9 

9 12 

3 9 


1 4 
3 8 

5 11 
9 14 

6 12 


1 5 

3 11 

5 13 

10 11 


1 6 
4 7 

6 8 
10 13 


2 7 

4 8 

6 10 
11 12 


2 8 
4 10 

6 14 

11 14 


The valnes for this design are: 


7l(t/jy,0) - 68.2336883 
L)(V?^,0) = 1,627,763,046 
A{'ijjT,5) - 1054.7827063 

5) = 2, 529, 608, 088, 061, 601, 720,121, 676. 


2 9 

4 12 

7 14 

12 13 


The R-best RGD for ?/ > 0 and D-best RGD for ?/ > 0 is #12 in the appendix; its valnes 
are 


R(#y,0) - 68.2334019 
L)(#^,0) = 1,627,920,000 


R(#y,5) - 1054.7827069 


Di'ipT, 5) = 2, 529, 608, 091, 727, 840, 200, 000, 000. 


Remark 2. Let denote the regnlar complete mnltipartite graph, that is a graph 

whose vertex set can be partitioned into gronps of size m snch that any pair of vertices is 
joined by an edge if and only if they are in different gronps. An RGD with a mnltipartite 
concnrrence graph and block size 2 is a gronp divisible design. 


The RGDs corresponding to the adjacency matrices of A" 2 , 2 , 2 , A" 2 , 2 , 2 , 2 , 772 , 2 , 2 , 2 , 2 , 772 , 2 , 2 , 2 , 2 , 2 , 
773 , 3,3 and 773 , 3 , 3,3 stay best for all x > 0 . 


Ghena fll98ll) proved that regnlar complete bipartite graphs give rise to the concnrrence 


graphs of the nniqne A- and D-optimal designs for all x > 0 (not necessarily only among 
RGDs). He extended his resnlt to complete regnlar mnltipartite graphs among simple 
graphs. 
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Remark 3. We found support fo r the f ollowing conjecture. 

Conjecture 1 ( Bagchi &: Chenj (1993)). The A- and D-best RGBs with block size 2 
have an adjacency graph that follows the pattern below. 


Take the complete bipartite graph with parts of size v/2 — 1 and v/2 + 1 and add on the 
larger part the edges of a circuit on n/2 + 1 vertices, for example for n = 10 and 5 = 6: 



These graphs are A- 
For n = 14 we could 


and H-best for n < 14 which was found bv iJohn fc Mitchell! (1l977l) . 
verify the A- and H-best for this graph, too. 
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A Optimal regular graph designs 


#1 u = 7, A; = 2, r = 4 + 6|/ 
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w = 
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k = 
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#21 V = 20, A; = 2, r = 3 + 1% 
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17 
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14 
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11 
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19 

11 

18 

17 

20 

18 

20 


3 8 4 9 

7 15 8 13 

12 16 13 20 
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